Disorder in Diluted Spin Systems 
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The influence of substitutional disorder on the magnetic properties of diluted Heisenberg spin 
systems is studied with regard to the magnetic stability of ferromagnetic diluted semiconductors 
(DMS). The equation of motion for the magnon Green's function within Tyablikov approximation 
is solved numerically for finite systems. The resulting spectral density is then used to estimate the 
magnetization and Curie temperature of an infinite system. This method is suitable for any form 
of a ferromagnetic exchange interaction. Besides different lattices and spin magnitude S, exchange 
interactions of different range are examined. The results show that, for short-range interaction, no 
magnetic order exists below the critical percolation concentration, whereas a linear dependence of the 
Curie temperature on the concentration of spins is found for ferromagnetic long-range interaction. 

PACS numbers: 75.10.Nr, 75.10.Jm, 75.50.Pp, 02.60.-x 



I. INTRODUCTION 

The discovery of ferromagnetism in (III,Mn)-V di- 
luted magnetic semiconductorsiiSi^ (DMS) has at- 
tracted considerable attention among scientists during 
the past years i^iSiSiSiiS The great interest in these ma- 
terials is motivated by the idea of using their spin 
degrees of freedom in conjunction with their elec- 
tronic degrees to build new electronic devices ranging 
from fast nonvolatile memories to quantum comput- 
ers (spintronics)iii*iSiiiiiiiS-i6'^'^ However, progress has 
been impeded by the fact that most of the DMS stud- 
ied so far have a Curie temperature Tq below room 
temperature2iiiLi2iiS. 

For the development of ferromagnetic DMS with de- 
sired properties such as high Curie temperatures, the 
theoretical understanding of the magnetism in these ma- 
terials plays an important role. In DMS, a small frac- 
tion of the non-magnetic host-semiconductor ions is sub- 
stituted by ions, which carry a localized magnetic mo- 
ment (spins). These magnetic ions are mainly ran- 
domly distributed over the lattice sites of one of the 
host-semiconductor species. This positional disorder 
breaks the translational symmetry of the crystal and thus 
greatly complicates the theoretical description of the ma- 
terial. 

One of the first to consider the substitutional 
disorder present in diluted spin systems was Brout 
in the late 1950si^ Since then, various attempts 
have been made to tackle the problem of includ- 
ing the disorder into the theoretical description 
of ferromagnetism in diluted spin systems. Many 

^j-.p^..„^„^„ j21.22.23.24.25.26.27.28.29.3 0.31 .32.33,34.35.36.37 

combine mean-field or spin-wave theories with 
configurational-average methods, which not only 
rely on uncontrolled approximations with respect 
to the disorder, but also encounter difficulties when 
applied to spin systems^^. Some use Monte-Carlo 
simulationa^^i-^'^i^^i^^'^'^i'^^'^^i^^i^^ which make the sim- 
plification of classical spins instead of quantum spins. 
Other approaches are based on percolation theorji^2i^2iS 



or replica methods^^'^^, but these treat the magnetism 
itself only on a mean-field level. However, one must 
properly take into account the positional disorder of 
the spins and their quantum nature to obtain reliable 
predictions about the magnetic properties of diluted 
spin systems 
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In this article, we will present an approach, which is 
based on the Tyablikov approximation^ to the Heisen- 
berg model and uses numerical studies of finite systems 
to estimate spontaneous magnetizations and Curie tem- 
peratures of the corresponding systems in the thermody- 
namic limit. It should be mentioned here that the idea 
of using numerical studies of finite spin systems in spin- 
wave-type approximation is not new^SiSiS According 
to our knowledge, however, the obtained spectral densi- 
ties have never been used before to calculate spontaneous 
magnetization and Curie temperatures. The presented 
approach is suitable for any form of ferromagnetic ex- 
change interaction. Within this approach, we are able 
to treat the quantum fluctuations of the spins within 
random-phase approximation. It means that our ap- 
proach goes beyond the classical-spin approximation and 
the mean-field theory, which notoriously overestimates 
ferromagnetic stability. Furthermore, we are able to 
treat the positional disorder of the spins numerically ex- 
act in contrast to approaches employing configurational- 
average methods. 

This article is organized as follows: In sectional we 
will generalize the Tyablikov approximation to spin sys- 
tems without translational symmetry. Section [llll is con- 
cerned with the numerical studies we performed on the 
basis of the generalized Tyablikov approximation. In sub- 
section nil Al we will explain the algorithm we used for 
our numerical studies. In the following subsections, we 
will discuss the results we obtained for systems with fer- 
romagnetic short-range, ferromagnetic long-range and os- 
cillating long-range exchange interactions. In scction llVl 
we will conclude the article with a summary and propos- 
als for future work. 
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II. THE MODEL 

To study the magnetic properties of a material with 
a concentration c e [0, 1] of atoms carrying a locahzed 
magnetic moment, we consider a lattice with a fraction 
c of the lattice sites occupied by a spin. We assume the 
dynamics of these spins to be described by the isotropic 
Heisenberg Hamiltonian 



N N 



(1) 



Here, i and j label the occupied lattice sites and 
Si = {Sf , Sf , Sf ) is the spin operator of the locahzed 
magnetic moment at lattice site i with lattice vector R^. 
The spins - whose total number is denoted by iV - in- 
teract via a Heisenberg exchange coupling with exchange 



parameters Jy obeying Jij — J(|Ri — Rj|) and Ju 



0. 



In addition, the Hamiltonian contains a Zeeman coupling 
of the spins to a homogeneous external magnetic field 
B = (0,0,5). 

The Tyablikov approximation™ has been developed for 
and very successfully applied to concentrated spin sys- 
tems with translational symmetry. In the following, we 
will generalize the Tyablikov approximation to derive ex- 
pressions for the magnetization and the Curie tempera- 
ture for diluted spin systems, where we can not presup- 
pose translational symmetry. Thereby, we will keep the 
form of the appearing expressions as close as possible to 
those known from the usual Tyablikov approximation. 
However, the absence of translational symmetry compli- 
cates the evaluation of the derived expressions, as will be 
seen later. 

We will first consider spins with S — 1/2. In this case, 
the magnetic properties of the system can be studied by 
use of the retarded magnon Green's function 



ret 
J // E 



(2) 



containing the step operators S.^ = Sf ± iSf. The equa- 
tion of motion for Gy (E) reads 

{E - gjti^B) G,j{E) = 2h'5,j (Sf) + J.^x 

m 

^ {^^t^mT^j )) B * ~ {{Sm^ii^j )) B *) • (3) 

Making the Tyablikov approximation, which consists 
in decoupling the higher Green's function on the 
rhs. of (Q, and assuming a uniform magnetization 
(S,;) = (0,0, (S*^)), one obtains after rearrangement: 



E 



E- gjUBB 



J it 



hSi 



. (4) 

In the following, it is convenient to interpret this equation 
as a matrix equation: 



1 - H G{E) = hi 



(5) 



Here, the Hamilton matrix H, which is defined by its 
matrix elements Hij = Sij X]^)!=i "^in ~ Jij, has been in- 
troduced. 

In order to solve equation ^ for G, one has to diago- 
nalize H. In the concentrated case, this can be done by a 
Fourier transformation because of the translational sym- 
metry of the system. For a diluted system, translational 
symmetry is absent, and hence, no general method exists 
for this task. However, for a finite system the Hamil- 
ton matrix can be diagonalized numerically. Since H is 
Hcrmitian, this diagonalization yields N real eigenvalues 
Er with eigenvectors, which can be chosen to form an 
orthonormal basis. For notational convenience, we will 
use indices r and r' to refer to matrix elements in this 
'H-diagonal' basis, whereas we retain i and j to refer to 
matrix elements in the 'site-diagonal' basis, where matrix 
indices are identified with site labels. 

After going to the H-diagonal basis, one obtains for 
the magnon Green's function: 



Grr'{E) — 



E - gj^BB - 2h (S^) Er + iO+ 



(6) 



From this expression, one can read off the magnon spec- 
tral density Srr'{E) = (G„/(£;)) /n: 



Srr'{E) 



-3(G„K^)) 



= 2h'^ (5") S{E- gjfiBB - 2h {S') Er) 6r 



(7) 



Thus, the energies Er = gjfiBB — 2h (S^) Er are the ex- 
citation energies of the elementary excitations of the sys- 
tems, which we call magnons. 

To obtain the local magnon spectral density Sui^E), 
one has to transform Srr'{E) back to the site-diagonal 
basis. However, this may yield different SuiE) for dif- 
ferent lattice sites i. This contradicts the assumption of 
a uniform magnetization, since Sii{E) and (Sf) are con- 
nected by the spectral theorem: 

(5f ) = hS^]:{SrS^) = hS-^ 1^ ^^dE , (8) 

where (3—1 /keT with ke as the Boltzmann constant and 
T as the temperature. One can circumvent the problem 
by abandoning the simplification of a uniform magneti- 
zation at the cost of having to solve a system of N non- 
algebraic equations in an iterative procedure, where one 
has to diagonalize an N x N matrix similar to H in each 
step. Here, a numerically less demanding approach is 
taken, which we believe to capture the important physics 
of the system: We consider (S^) an approximation for 
(5,f ) and use the site-averaged spectral density 

N N 

Sh{E) = ^ E ^"(^) = ]^ E ^-(^) (9) 
as an approximation for 5'^^ {E) and then use the spectral 
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theorem: 



Su{E) 



dE 



"S-^J'^.^E. ,10) 



In short, we have to solve: 



SujE) 



dE 



(11) 



After inserting ^ and doing some rearrangements, this 
yields: 



hs 



-■ (12) 



From this equation, one obtains an expression for 
the Curie temperature by linearization with respect to 



Tn = 



2kR I TV ^ Er 



(13) 



Note that for translationally invariant systems, the in- 
dices r can be identified with the wave vectors of the 
first Brillouin zone, and thus the expressions H12|l and 
(|13|l reduce to the standard Tyablikov formulae for con- 
centrated spin systems. 

For general S'-values, an analogous derivation using a 
set of Green's functions^^'^^'^^'^^ instead of the Green's 
function |(2J) alone leads to the implicit equation 



{SI = 



(1 + 5' + $) $^^+^ + (5* - (1 + $) 
(1 + $)^'^+^ - $2S+1 



2S+1 



(14) 

for the magnetization. Here, the average magnon number 



N ^ 



1 



TV ^ f.l3{gjiisB+2h{S')Er.) _ 1 



(15) 



has been used as a convenient abbreviation. The Curie 
temperature is then given by 



(16) 



for general S'-values. 

In the following we will only consider systems in an 
infinitesimal external magnetic field B = 0"*" just to break 
rotational invariance and therefore, usually drop terms 
including the i?-field explicitly. 



III. NUMERICAL STUDIES OF FINITE 
SYSTEMS 

A. Algorithm 

The equations (|12|l to Hlt)|) are of no practical use for 
the calculations of the magnetic properties of the sys- 
tem unless one can compute the eigenvalues Er of the 
Hamilton matrix H. As already mentioned above, the 
Fourier transformation, which is the weapon of choice 
for concentrated spin systems, does not work for diluted 
systems, because translational symmetry is absent. How- 
ever, for a finite system, we can diagonalize the Hamilton 
matrix numerically to obtain its eigenvalues. Compared 
to configurational average methods, the numerical diag- 
onalization has the advantage that no approximations 
with respect to the disorder are needed. Furthermore, 
the spins can be treated quantum mechanically and the 
simplification to classical spins as in usual Monte Carlo 
simulations is not necessary. 

The algorithm we used in our numerical studies of fi- 
nite systems is as follows: First of all, the system parame- 
ters must be fixed. These are the concentration c of spins 
on the lattice, the exchange parameter J (R) as a func- 
tion of the distance R between two spins and the size and 
geometry of the system. The next step is to decide by the 
help of a pseudo-random-number generator which lattice 
sites of the system are occupied by a spin. Now that 
the positions of all occupied lattice sites in the system 
are known, the exchange parameters Jy = J (|Ri — R.j|) 
between occupied lattice sites i and j can be calculated, 
and hence, the matrix elements Hij . Now, the matrix 
H can be diagonalized numerically using standard diag- 
onalization routines for large matrices. 

For the shape of the simulated samples, we chose a 
cuboid lattice section, which we closed to a torus to 
eliminate surface effects. To avoid multiple- and self- 
interaction due to the periodic boundary conditions, we 
cut off any long-range interaction at a distance half the 
system size. We chose the size of the samples such that 
for each concentration we obtained approximately the 
maximal number of spins in the system that we could 
handle (ca. 14 000 spins on an AMD Athlon PC, 2GHz 
CPU, 1 GB RAM). For simplicity, we only considered 
the case that each lattice site is occupied by a spin with 
probability c independent of the occupation of neighbor- 
ing sites. However, if one is interested in clustering ef- 
fects, a different routine to distribute the spins on the 
lattice could be used. 

To decrease statistical errors, we averaged the eigen- 
value distribution, i.e. Sii{E), over several simulated 
systems with equal system parameters, but different 
pseudo-random distribution of the spins. Test calcu- 
lations showed that this averaging seems hardly neces- 
sary for the large samples we used, which is due to the 
self-averaging property of the site-averaged local magnon 
spectral density Sii{E)^. As a trick to smoothcn the 
eigenvalue spectra, we sometimes used slightly differ- 
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ent sample geometries for different simulation runs and 
weighted the obtained eigenvalue distributions with their 
number of lattice sites^. Further weighting was not nec- 
essary, since for a fixed system size the probability of 
generating a certain distribution of spins with the used 
algorithm coincides with the statistical weight of this dis- 
tribution among those having an equal number of lattice 



sites 



64 



The disadvantage of the numerical diagonalization is 
that only finite systems can be treated, but there is no 
spontaneous magnetization in the isotropic Heisenberg 
model for any finite system^^*S& The Tyablikov approxi- 
mation preserves this property, which is seen in our cal- 
culations by the fact that, due to its structure, H has 
at least one zero-eigenvalue, which has a finite spec- 
tral weight for finite systems. More precisely, the num- 
ber iVzoios of zero-eigenvalues of H equals the number 
-^clusters of clustcrs of spius Connected by a non-zero ex- 
change interaction for systems with only ferromagnetic 
interaction and Nzcms > -^clusters for systems contain- 
ing antiferromagnetic interaction.^^ Hence, (|12|) and 113() 
always yield zero for the spontaneous magnetization and 
the Curie temperature, if the computed eigenvalues of the 
simulated samples are used directly as input. In order to 
obtain a finite spontaneous magnetization and Curie tem- 
perature, one has to go to infinite systems. However, the 
eigenvalue distribution of a sufficiently large finite sys- 
tem represents a good approximation to the eigenvalue 
distribution of the corresponding infinite system except 
for energies close to zero, where the discreteness of the 
eigenvalue spectrum of the finite system comes into play. 
We solved this problem by shifting the zero-eigenvalues 
of H slightly to higher energies, which will be explained 
in more detail in the following subsection. 



B. Short-Range Interaction 

As first application we investigate systems of spins on 
a simple cubic (sc) lattice with nearest-neighbor interac- 
tion only: 



JiR) 



jO for i? = a, 
otherwise. 



(17) 



Here, a denotes the nearest neighbor distance and J*' the 
nearest neighbor interaction strength. 

In figure m we plotted the spectral density for differ- 
ent concentrations c of spins. The spectral density is 
zero for all i? < for all concentrations, which is due 
to the fact that H is positive semi-definite for systems 
with purely ferromagnetic interactions. The calculations 
for the concentrated system show the well-known sym- 
metric shape of a simple cubic density of states. Dilu- 
tion increases the spectral density for lower energies at 
the cost of the spectral density at higher energies. For 
concentration above c ~ 0.5, the spectral density stays 




FIG. 1: Site-averaged local magnon spectral density Sii{E) 
for nearest-neighbor interaction on an sc lattice for various 
concentrations c of spins (shown as histogram of the eigen- 
value distribution of the Hamilton matrix, each curve is an 
average of 5 to 10 configurations of 8 000 to 14 000 spins) 



smooth, whereas for lower concentrations, peaks at cer- 
tain energies {E/2H {S^) J° = 0, 1, 2, 2 - v^, . . .) become 
apparent, which constitute the whole spectral density for 
c < 0.3. The smooth part of the spectral density can be 
attributed to one large percolating cluster of spins, which 
has a quasi-continuous spectrum. This percolating clus- 
ter only exists above the critical percolation concentra- 
tion cp, which is Cp ~ 0.3 for the simple cubic lattice with 
nearest-neighbor interactioniSSiSSi The peaks stem from 
smaller clusters, which contain a considerable fraction of 
the spins only for concentrations around and below cp 
and which have a discrete spectrum. 

To obtain non-zero Curie temperatures from the ob- 
tained H-eigenvalues of finite systems, we used a small 
energy 6 to shift the eigenvalues close to zero slightly to 
higher energies. To be more specific, we used 




for Er < S, 
otherwise 



Vr e {1,...,7V} (18) 



for practical reasons. The effect of this shift 6 on the 
Curie temperature Tq is illustrated in figure |21 Let us 
first discuss the curves for c = 1, where the numeri- 
cally obtained data for finite systems can be compared 
with the analytical results for the infinite system. For 
N = oo, the Curie temperature has a finite value at 
S ~ Q and is quite insensitive to small shifts (5 ^ J'^. 
For large enough but finite N, the calculated Curie tem- 
perature rapidly rises from Tc = at 5 = to a value 
comparable to the infinite system at (5 w /N . In this 
steep region, Tc{5) is dominated by the shifted zero- 
eigenvalue, which has a spectral weight Wzeros = N~^, 
and Tc{6) w ^fi^S {S + 1) N6. For S > J°/N, however, 
the shifted zero-eigenvalue becomes less important, and 
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the curves of the finite system and the infinite system are 
almost identical. Hence, the Curie temperature for the 
finite system with c = 1, large N and /N ^ (5 <C J" 
approximates the Curie temperature for the infinite sys- 
tem with c — \ very well. Furthermore, for large N and 
J'^/N ^ S ^ J^, the results are quite insensitive to the 
actual choice of S. 

For c < 1, the situation is more complicated. Firstly, 
there are no exact analytical results for Tc we could com- 
pare our numerical results with. Furthermore, for an 
exchange interaction with finite range, there is always 
a finite fraction of spins in smaller clusters. Hence, al- 
ways a finite fraction Wzcms = Nzeros/N of eigenvalues is 
zero. Therefore, Tc{S) maintains a finite slope oc Wzeros 
for 6 ^ WzcTosJ^ even for very large N as is seen in 
figur^b). This problem could be avoided if only the 
percolating cluster, which is the only cluster that can 
support long-range order in the system, were considered 
in the calculation of Tc. However, identifying and remov- 
ing all the eigenvalues due to smaller clusters is numeri- 
cally quite demanding. Moreover, for concentrations well 
above the critical percolation concentration, the number 
of smaller clusters is very smallS and thus Wzeros <C 1. 
Hence, the behavior of Tc{S) is dominated by the eigen- 
values of the percolating cluster for WzemsJ'^ S and 
c > cp and is not very sensitive to changes in S for 5 ^ j" 
as seen in figure |2Jb). 

For c < Cp most or all eigenvalues are due to smaller 
clusters and therefore, the applicability of the presented 
method may be questioned. However, since a large frac- 
tion of the eigenvalues is zero, Tc{S) « ^h'^S{S + 1)6, 
which simply means that Tc ~ for small S. This 
agrees with percolation theorj^, from which follows that 
Tc = for c < Cp, since there can not exist any long- 
range ordering in the system. Therefore, we believe that 
the Curie temperature calculated using the eigenvalues 
of a large finite system and a shift 6 chosen such that 
J°/N <C S <C J° gives a good estimate for the Curie 
temperature of the corresponding infinite system for fi- 
nite dilution, too. In the following, all given values for 
the Curie temperature and spontaneous magnetization 
refer to S = 0.01J°, which corresponds for systems with 
N « 14 000 to a Curie temperature Tc{S) just outside 
the steep region at small S. 

The Curie temperature Tq we calculated for various 
concentrations is shown in figure |3 Below c — 0.3, 
the Curie temperature is practically zero in accordance 
with percolation theory. Above c = 0.3, the Curie 
temperature increases with increasing concentration and 
reaches its maximum value keTc — 2.615(5' + l)h^j'^ at 
c = 1 in good agreement with the analytical result of 
keTc = 2MS{S+1)H^J°. 

The question of the effect of the substitutional disor- 
der at finite dilution on the Curie temperature is con- 
nected to the question, how one defines the ordered sys- 
tem for finite dilution. Placing the remaining spins of 
the diluted system on an artificial lattice with an appro- 
priate lattice constant to regain translational invariance 
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FIG. 2: Curie temperature Tc as function of the energy shift 
S used for the low-lying eigenvalues of H for different number 
A'' of spins in the system for (a) c = 1 and (b) c — 0.7 
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FIG. 3: Curie temperature Tc as function of the concentra- 
tion c for nearest-neighbor interaction on an sc lattice 



would give a Curie temperature merely depending on 
how one alters the nearest-neighbor exchange interaction 
strength with increasing lattice constant. The virtual- 
crystal approximation, which is widely used to incorpo- 
rate dilution effects because of its simplicity, would give 
Tc(c) — crc(I), which is certainly wrong for low c in the 
case of an exchange interaction with finite range. 

The reduced magnetization as a function of tempera- 
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FIG. 4: Reduced magnetization {S') / (S^) (T = 0) as func- 
tion of reduced temperature T/Tc for nearest-neighbor in- 
teraction on an sc lattice for various concentrations c for (a) 



taken into account for the spin-spin interaction in the cal- 
culations), a is again the nearest-neighbor distance and 
J° the nearest- neighbor interaction strength. For n = 1, 
we recover the case of pure nearest-neighbor interaction, 
which we discussed in the last subsection. For n — > oo, 
we have a long-range interaction with algebraic fall-off. 
In this case, there is always only one cluster even for fi- 
nite dilution, and we do not encounter complications due 
to smaller clusters in the calculation of the Curie tem- 
perature. 

In figure El we plotted the Curie temperature as 
function of the concentration c for different ranges n of 
the exchange interaction. We performed calculations for 
the sc and fee lattice structure, but no qualitative differ- 
ences could be found between them. As expected, a de- 
creasing value of the threshold concentration, above one 
finds a finite Curie temperature, is found with increas- 
ing range n of the interaction. Moreover, the observed 
values of the threshold concentration agree quite well 
with the corresponding critical percolation concentra- 
tion, which is a lattice-structure- and interaction-range- 
dependent quantity. ^-'-^'^^ Furthermore, the results show 
that the exchange interaction must have a range beyond 
the fourth shell to obtain a finite Curie temperature for 
concentration values of c ~ 0.05, which are typical of the 
recently studied DMS. This is consistent with the values 
n > 5 for the fee and n > 6 for the sc lattice obtained by 
the formula 



S = 1/2 and (b) S = 5/2 



n> ( —f 



2/3 



(20) 



ture is shown in figure 0] for several concentrations. For 
lower concentrations, the magnetization curves show a 
more linear decrease with temperature. This effect is less 
pronounced for higher spins. Any anomalous curvature 
can not be found. 



C. Ferromagnetic Long- Range Interaction 

The experimental fact of Curie temperatures of more 
than 100 K in DMS such as Ga-^ .^.Mn^As with concen- 
trations as low as ~ 0.05 of magnetic atom^iiSSiSi implies 
that the exchange interaction between the localized mo- 
ments is not only quite strong, but must be very long- 
ranged in these materials. To study the influence of the 
range of a ferromagnetic exchange interaction on the sta- 
bility of the ferromagnetic phase with respect to dilution 
we choose 



J(i?) 




for a < R < ^/na, 
otherwise. 



(19) 



Here, n is the number of the outermost shell with non- 
zero exchange interaction (or the outermost shell that is 



with /|e'om = 2.4 and f^Z^^ = 2.4/^2, which is based on 
percolation arguments 4^ 

For n — oo, the Curie temperature is essentially linear 
in the concentration for both lattice structures. Hence, a 
small, but finite value is found for the Curie temperature 
even for very low concentrations. Thus, it seems that 
the applicability of the virtual-crystal approximation is 
restored in the limit n — > oo. However, calculation for 
other forms of ferromagnetic long-range exchange inter- 
actions show that the perfect linear dependence T^^ oc c is 
a peculiarity of the exchange interaction l|19|l . although 
the proportionality seems to be roughly fulfilled for all 
ferromagnetic long-range exchange interactions. Note 
that placing the remaining spins of the diluted system on 
an artificial lattice with an appropriate lattice constant 
a(c) cx c^^^ and changing the exchange parameter accord- 
ing to eq. H19(l gives a Curie temperature T^-, cx c^^^ in 
contrast to the observed cx c. 

We also calculated magnetization curves for different 
concentrations and exchange-interaction ranges for the 
sc and fee lattice structure. However, since any unusual 
features of the magnetization curves have not been ob- 
served, we will refrain from an extensive discussion of 
these. 



7 




0.2 0.4 0.6 0.8 1 

c 




0.2 0.4 0.6 0.8 1 



c 

FIG. 5: Reduced Curie temperature Tcic) jTciV) as function 
of the concentration c for exchange interactions with different 
range n for (a) sc and (b) fee lattice structure 



D. Oscillating Long- Range Interaction 

It is believed that the ferromagnetism in DMS is 
carrier-induced, with holes in the valence band of the host 
semiconductor mediating the interaction between the lo- 
calized moments of the impurity atoms. One feature of 
such an indirect exchange mechanism is that it may give 
rise to a long-ranged effective exchange interaction J(-R) 
which is non-monotonous as function of the spin-spin dis- 
tance i? and even may be ferromagnetic for some, but 
antiferromagnetic for other values of R. To study such 
an oscillatory exchange interaction, we choose 

J(i?) = P°(?)^'^°<"(^"l)) (21) 
1 otherwise. 

In figure IHl we plotted the site-averaged local magnon 
spectral density S^^{E) for different concentrations. Al- 
though the exchange interaction l|21|l contains antifer- 
romagnetic interactions, the concentrated spin system 
has a saturated ferromagnetic ground state and a low- 
temperature ferromagnetic phase with a Curie tempera- 
ture of kgT(^ « 1.6fi.^J°. However, as soon as the con- 
centration is decreased from c = 1, the spectral density 
becomes non-vanishing for £' < 0. This effect is still small 



c= 1 




EI2h{S~)f 



FIG. 6: Site-averaged local magnon spectral density Sii{E) 
for oscillating long-range exchange interaction on an sc lattice 
for various concentrations c of spins (shown as histogram of 
the eigenvalue distribution of the Hamilton matrix, each curve 
is an average of 5 to 10 configurations of 8 000 to 14 000 spins) 



for high concentrations, and one could try to use the com- 
puted eigenvalues to calculate Curie temperatures by us- 
ing a small shift of the lowest eigenvalues, but already at 
intermediate concentrations the negative H-eigenvalues 
substantially contribute to S^^ and certainly can not be 
ignored any more. This means that we can not compute 
any spontaneous magnetization or finite Curie tempera- 
ture for such systems with our method. 

The results observed for the exchange interaction 1)2111 
can be generalized. Eigenvalues i?,, < of the Hamilton 
matrix H occur for all systems with c < 1 and antiferro- 
magnetic components in the exchange interaction. These 
negative eigenvalues stem from the fact that at finite di- 
lution there are always some spins in the system which 
predominantly couple antiferromagnetically to the other 
spins of the system. As a result, the saturated ferromag- 
netic state is not the ground state of the diluted system 
and hence, is thermodynamically unstable against cre- 
ation of magnons. However, the Tyablikov approxima- 
tion with the assumption of a uniform magnetization is a 
reasonable approximation to the Heisenberg model only 
for a saturated ferromagnetic ground state, and therefore, 
the presented method fails for oscillatory exchange inter- 
actions. Note that these considerations rule out a highly 
oscillating exchange interaction for materials which are 
ferromagnetic despite a low concentration of magnetic 
atoms. For materials with highly oscillating exchange in- 
teraction, one rather expects a spin-glass-like phase at 
low concentrations instead of a ferromagnetic one^S 
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IV. SUMMARY 

The aim of this article was to study the magnetic prop- 
erties of diluted spin systems. Our first step was to gener- 
alize the Tyablikov approximation to spin systems with- 
out translational symmetry to account for the substitu- 
tional disorder present in such systems. The resulting 
equation of motion for the magnon Green's function was 
then solved numerically for finite spin systems. The ob- 
tained spectral densities were then used to estimate the 
Curie temperature and the spontaneous magnetization of 
the corresponding infinite system. 

Compared with other 

^,. 21. 22. 23. 24. 25. 26. 27.28. 29. 30. 31. 32. 33. 34. 35. 36. 37. 39.40. 41.42. 43 

our approach has the advantage that the spins are treated 
quantum-mechanically within the Tyablikov approxima- 
tion, which goes beyond mean-field and classical-spin 
approximation, and - apart from the simplification of a 
uniform magnetization - no approximations are needed 
with respect to the substitutional disorder. Furthermore, 
the numerical effort as well as the analytical is fairly low 
in our approach. 

Our calculations for short-range exchange interactions 
show no magnetic order below the critical percolation 
concentration in accordance with percolation theory. For 
ferromagnetic long-range exchange interactions, the cal- 
culation show a linear dependence of the Curie tempera- 
ture on the concentration of spins. Hence, a finite Curie 
temperature is found even for very low concentrations. 
For systems with oscillating long-range exchange inter- 
action, the presented method fails to give answers due 
to the lack of a ferromagnetic ground state at finite dilu- 
tion. For ferromagnetic materials with concentrations of 
magnetic atoms as low as c « 0.05, these results imply 
that the effective exchange interaction must be very long- 



ranged, but can not oscillate strongly with the inter-spin 
distance. 

In this article, we have shown that our approach gives 
quite reasonable results for a variety of model systems 
with ferromagnetic ground state. The next step should 
be an application to real substances. To attack the prob- 
lem of ferromagnetic DMS, for example, the presented 
method may be combined with a theory for the elec- 
tronic degrees of freedom, which provides the values for 
the effective exchange interaction between the localized 
moments. Furthermore, the method can be extended to 
account for a site-dependent magnetization of the local- 
ized moments. However, this will highly increase the nu- 
j|t*'*'*'i*>l'^'^^^^'5SltS& Clustering and other forms of chemical 
ordering may be easily included into the model to make 
it more realistic for materials, in which these play a role. 
Moreover, the method could then be used to analyze of 
the effects of chemical short-range order on the magnetic 
properties. This knowledge could be very useful for the 
development of new materials, in which one can influ- 
ence the distribution of the magnetic atoms by choosing 
certain preparation parameters and hence, obtain certain 
desired magnetic properties. Furthermore, it should be 
investigated whether the presented method, in particular 
the treatment of the zero-eigenvalues, could be put on a 
firm theoretical basis. 
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